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Abstract. After a short introduction to the nonmesonic weak decay (NMWD) A^V nN of A- 
hypernuclei we discuss the long-standing puzzle on the ratio Fn/Fp, and some recent experimental 
evidences that signalized towards its final solution. Two versions of the Independent-Particle-Shell- 
Model (IPSM) are employed to account for the nuclear structure of the final residual nuclei. They 
are: (a) IPSM-a, where no correlation, except for the Pauli principle, is taken into account, and (b) 
IPSM-b, where the highly excited hole states are considered to be quasi-stationary and are described 
by Breit-Wigner distributions, whose widths are estimated from the experimental data. We evaluate 
the coincidence spectra in ^He, ^He, ]^C, ]^0, and ^^Si, as a function of the sum of kinetic energies 
E„N = E„ +En forN — n,p. The recent Brookhaven National Laboratory experiment E788 on ^He, 
is interpreted within the IPSM . We found that the shapes of all the spectra are basically tailored 
by the kinematics of the corresponding phase space, depending very weakly on the dynamics, 
which is gauged here by the one-meson-exchange-potential. In spite of the straightforwardness of 
the approach a good agreement with data is achieved. This might be an indication that the final- 
state-interactions and the two-nucleon induced processes are not very important in the decay of 
this hypernucleus. We have also found that the k + K exchange potential with soft vertex-form- 
factor cutoffs (A^ w 0.7 GeV, Ak « 0.9 GeV), is able to account simultaneously for the available 
experimental data related to Fp and F„ for ^H, ^He, and ^He. 
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INTRODUCTION 

The nonmesonic weak decay (NMWD) of A hypernuclei, AiV nN {N = p,n), takes 
place only within nuclear environment. Without producing any additional on-shell par- 
ticle (as does the mesonic weak decay A — )■ tzN) the mass is changed by 176 MeV, 
and the strangeness by AS = 1, which implies that we are witnessing the most drastic 
metamorphosis of an elementary particle within the nucleus. As such, the hypernuclei 
can be considered as a powerful "laboratory" for unique investigations of baryon-baryon 
strangeness- changing weak interactions, and the NMWD could play an important role 
in the stability of rotating neutron stars with respect to gravitational wave emission [1,2]. 

Same as the free A hyperon, they are mostly produced via the strong interactions, 
i.e., in the reaction processes Tt^n — )• AX+, K n — )• K^A and K p — )■ Tt^A, by making 



use of the pion (n) and kaon (K) beams. They also basically decay through the weak 
interactions, as the free A does. Yet, as it is well known and explained below, there are 
some very important differences in the corresponding decaying modes. 




FIGURE 1. On the left and right panels are represented, respectively, the standard nuclear chart of 
nuclides, and the hypernuclear chart of isotopes. The NMWD corresponds to a jump from the surface 
S=l in the three-dimensional chart (N,Z,S) to the standard (N,Z) chart. 

Not less important is the fact that with the incorporation of strangeness, the radioac- 
tivity domain is extended to three dimensions (N,Z,S), as illustrated in Fig. 1. The 
best studied systems are nuclei containing a A-hyperon, which because of the additional 
binding are even richer in elements than the ordinary (N, Z) domain. (For instance, while 
the one-neutron separation energy in ^^C is 1.01 MeV, it is 1.63 MeV in ^^C [3].) 
This attribute of hypernuclei has motivated a recent proposal to produce neutron rich 
A-hypernuclei at J-PARC, including ^He [4]. The shrinkage of the ^'^C nucleus by the 
addition of an A-hyperon to build up the ^^C hypemucleus is illustrated in Fig. 2. 




FIGURE 2. Pictorial representation of shrinkage of C nucleus by the addition of an A-hyperon to 
build up the ^^C hypernucleus. 

The free A hyperon decays (as represented schematically by the first graph in Fig. 3) 
nearly 100 % of the time by the A — )■ Nn weak-mesonic mode 

j p + 7t (64.1%) 
^\ n + 7t^ (35.7%), 



with the total transition rate + F^q = F*^ = 2.50 ■ 10 ^ eV (which corresponds to the 

lifetime = 2.63 • 10~^^ sec). For the decay at rest the energy-momentum conservation 
implies 

MK = MM+^^ + \Jpl + ml; pN = Pn- 

Therefore the energy released is 

Qo = Ma -Mat -m;r ^ 37 MeV, 
and the kinetic energies and momenta in the final state are: 

2Ma 

- eo-7]v = 32MeV, 



PN = Pn=^ {Tn+MnY -M2 ^ 100 MeV/c. 

During the decay the isospin is changed by AT = 1/2 and 3/2 and its projection by 
AMj- = — 1/2. However, as the above experimental data can be accounted for fairly well 
by neglecting the AT = 3/2 component, one end up with AT = 1/2 rule, which leads to 
the estimate /F^o = 2, while the experimental result is 64. 1 /35.7 = 1 .80. 

Assuming the AT = 1/2 rule, the phenomelogical weak Hamiltonian for the process 
depicted in Fig. 31 can be expressed as: 

Hw = -/GFm|v7^(A;r+5;r75)0j,-TV/A(^ J ^ , (1) 

where Gptn^ = 2.21 x 10~^ is the weak coupling constant. The empirical constants 
Ayi = 1.05 and Byt = —7.15, adjusted to the observables of the free A decay, determine 
the strengths of parity violating and parity conserving amplitudes, respectively. The 
nucleon, A and pion fields are given by i///v and i//a and ^j^, respectively, while the 

isospin spurion ^ ^ ^ is included in order to enforce the empirical Ar = 1/2 rule. 

The free A hyperon weak decay is radically modified in the nuclear environment 
because the nucleon and the hyperon now move, respectively, in the mean fields Uj^ 
and Ua, which come from the A^A^ and A^A interactions. Un and Ua are characterized by 
the single particle energies (s.p.e.) ej^ and £a, and we have to differentiate between: 

• Mesonic Weak Decay (MWD): The basic process is again represented by the first 
graph shown in Fig. 3, and described by the hamiltonian (1). Yet, the energy-momentum 
conservation is different: 

MA = MAf-eA + ei+^P^|p^ + y^2+^, (2) 

where A is the mass number, and are the s.p.e. of the loosely bound states above the 
Fermi energy e^. They are of the order of a few MeV, while Ea is the energy of the 
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FIGURE 3. Schematic representation of: a) mesonic (nonleptonic) decay N + n induced by the 
weak vertex Hw, and b) nonmesonic weak decay AA'^ nN, engendered by the product of a weak vertex 
Hw and a strong vertex Hg. 
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FIGURE 4. a) One-nucleon, and b) two-nucleon induced A-decay in nuclei. Strictly speaking the OME 
diagrams are only valid for nonstrange-mesons n,p,(0, and tj . For strange-mesons K, and K* the vertices 
W, and S should be exchanged, i.e., W ^ S, and S ^W. 



Osi/2 State and goes from -11.7 MeV for to -26.5 MeV for Pb [5]. Thus, the 
corresponding Q-values 

QM^MA-MN-mj, + eA-el, (3) 

are significantly smaller than Qq, particularly for medium and heavy nuclei. This small 
value of Qm makes, as illustrated in [6, Fig. 2], the MWD to be hindered due to the Pauli 
principle. In fact, the experimental decay rates + F^o = Tm = Fm(A Nk) are of 
the order of F^ only for nuclei with A < 4, and they rapidly fall as a function of nuclear 
mass. For instance, in ]^ C: F^o/F*^ = 0.217 ±0.084 and F^-/fO = 0.052+g;gg. (For a 
recent theoretical study of the MWD see Ref. [7].) 
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FIGURE 5. (a) direct, and (b) crossed uncorrelated, and (c) correlated two-pion-exchange diagrams for 
the weak AA^ nN transition amplitude. 



• Nonmesonic Weak Decay (NMWD): New nonmesonic decay channels AN — )■ nN he- 
come open inside the nucleus, where there are no pions in the final state; it is represented 
schematically by the second diagram in Fig. 3. The corresponding transition rates can 
be stimulated either by protons, Tp = r{Ap — t- np), or by neutrons, r„ = r{An — )■ nn). 
The energy-momentum conservation and the Q- value are, respectively: 

2 2 I |2 

and 

Qnm = Ma - Ma, + Ea + 4- (5) 

Since the mean energy of the bound single-particle states is 4 ~ ~30 MeV, the Q- 
value is Qnm ~ 120— 135 MeV, and this is basically the kinetic energy of the two 
particles that are ejected from the hypernucleus. Therefore, the NMWD possesses a 
large phase space in the continuum, as illustrated in [6, Fig. 3], and the momenta p„, 
and Pa, of two outgoing nucleons are relatively large (~ 420 MeV). Therefore, the non- 
mesonic mode is not blocked by the Pauli principle, and dominates over the mesonic 
mode for all but the ^•-shell hypemuclei. 

It is assumed very often that the hypemuclear NMWD AN — )■ nN is triggered via the 
exchange of a virtual meson, and the obvious candidate is the one-pion-exchange (OPE) 
mechanism, where the strong Hamiltonian 

Hs = igNNKWNjsTt-tWN: (6) 

with ^A'A'^r = 13.4, accompanies the weak Hamiltonian (1) [8]. 

Later on, the full one meson-exchange (OME) model has been introduced by Dubach 
et al. [10], as schematically represented by the first graph in Fig. 4. Also are considered 
frequently the two-nucleon induced NMWD, represented by the second diagram in in 
Fig. 4. Here, one or two bound nucleons are expelled to the continuum by the nuclear 
ground state correlations, and one of them, together with the hyperon A, exchanges one 
meson giving rise to three decaying nucleons, i.e., AnN — )■ nNN. The corresponding 



decay rate is denoted as Tj, and the total weak decay rate of a A-hypernucleus is then: 



— Fyv/ + Tp^M, (7) 

where: 

Tm = r^j- + r^o, TjyM = Ti + r2, Fi = r„ + r^. (8) 

The OME potential is sometimes complemented with the contributions of uncorrelated 
{2n) and correlated (2n/o) two-pion-exchange [9], which are illustrated in see Fig. 5. 



Tn/Tp PUZZLE 

Large experimental values of the ratio r,,/^ = Tn/Tp in ^He and ]^C, measured before 
the year 2003 [11-13], were a cumbersome puzzle for the theorists during almost 
two decades, as schematically represented in Fig. 6. In fact, following the pioneering 
investigations of Adams [8] several calculations have been done within OPE coupling 
scheme of the total NMWD rate, and the ratio Tn/p = Fn/Fp reproducing reasonably 
well the first one, but failing badly for the second observable, (see Refs. [14-18], and 
references therein). 
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FIGURE 6. r„/rp puzzle: The experimental values of the ratio r„/p are: a) ^He T„/p ~ 0.93 ±0.55 



+ 1.12 



(BNL) [11], 1.97±0.67 (KEK) [12], 0.45 ± 0.1 1 ± 0.03 (KEK-E462) [49], andb) j^-^C: r„/p 1.33_o 



(BNL) [11], 1.87±0.59+y:i]^ (KEK) [12], 1.17+^:^^+^:f^ (KEK) [13], 0.51 ±0.13±0.05 (KEK) [50, 51], 
while the theoretical estimates for r„/p vary between 0.09 and 0.70 for ^He, and between 0.08 and 0.50 



for i^c 



The deficiency of the OPE model was attributed to effects of short range physics, 
which should be quite important in view of the large momentum transfers involved. 
Although there have been some attempts to account for this fact by making use of quark 
models to compute the shortest range part of the transition potential [19-23], most of 
the theoretical work opted for the addition of other, heavier mesons in the exchange 
process [6, 10, 24-38]. None of these models gives fully satisfactory results. Inclusion 
of correlated two-pion exchange has not been completely successful either [39, 40]. Nor 



have the addition of uncorrelated two-pion exchange, two-nucleon induced transitions 
or medium effects, treated within the nonrelativistic [15, 41^5] or relativistic [46] 
propagator approaches, been of much help. 

Yet, several important experimental advances in NMWD have been made in recent 
years, which have allowed to establish more precise values of the neutron- and proton- 
induced transition rates r„ and Tp, solving in this way the long-standing puzzle of the 
branching ratio r„/p. They are: 1) the new high quality measurements of single-nucleon 
spectra Sn{E), as a function of one-nucleon energy En = E done in Refs. [51-54], and 2) 
the first measurements of the two-particle-coincidence spectra as a function of the sum 
of kinetic energies En+E^ = E, SnN{E), and of the opening angle 0„Af = 0> SnN{cos 6), 
done in Refs. [49, 50, 54-57]. 



TRANSITION RATES 

To derive the NMWD rate we start from the Fermi Golden Rule. For a hypernucleus (in 
its ground state with spin Jj and energy Ejj) to residual nuclei (in the several allowed 
states with spins Jp and energies Ea^jp) and two free nucleons nN (with total spin S and 
total kinetic energy EnN = En + Eff), reads 

Fn = 2n [ \{p„pNS;aNJF\V\Ji)\^5iAa^j,-ER-E„N) fj^" (9) 

SaNJp ^ ' ^ ' 

where for the sake of simplicity we have suppressed the magnetic quantum numbers. 
The NMWD dynamics, contained within the weak hypernuclear transition potential V , 
will be described by the OME model, whose most commonly used version includes 
the exchange of the full pseudoscalar {%,K,r\) and vector {p,(0,K*) meson octets 
(PSVE), with the weak coupling constants obtained from soft meson theorems and 
SU{6)w [10, 31]. The wave functions for the kets |p„p^5M57/^M/7) and \JiMj) are 
assumed to be antisymmetrized and normalized, and the two emitted nucleons n and 
are described by plane waves. Initial and final short range correlations are included 
phenomenologically at a simple Jastrow-like level, while the finite nucleon size effects 
at the interaction vertices are gauged by monopole form factors [31, 34]. Moreover, 

r IPk + PnI^ EnN + 2cOS OnNy/EnEN 

^^^2M(A-2)- 1^2 ' ^^^^ 

is the recoil energy of the residual nucleus, and 

Aa^j, =AM + Ejj- Ea^jp , with AM = Ma - M = 176 MeV, (11) 

is the liberated energy. 

It could be convenient to perform a transformation to the relative and cm. momenta 
(P = 5(Pn - Pa')' P = P« + Pn), coordinates (r = r„ - r^v, j{r„ + r^v)) and orbital 
angular momenta I and L, and to express the energy conservation as 



EnN + Er— Aa^jp =ep + ep— Aaj^jp = 0, 



(12) 



where 

p2 p2 A A 
en=— , Er= — z -, ep= = -En (13) 

are, respectively, the energies of the relative motion of the outgoing pair, of the recoil, 
and of the total cm. motion (including the recoil). 

Following the analytical developments done in Ref. [34], the transition rate can be 
expressed as a function of the cm. energy ep: 

16M^ f A — 2\^^^ I 

Tat = -^[—^) / depY, ^ep{^a^J,-ep)^a^J,{pP). (14) 

It is understood that the square root should be replaced by zero whenever its argument 
is negative. Here 



SklLTJ 



and 



^ {plPLXSJT; jAjNJ^AN) {Ji 1 1 (^Af'^L) 1 1 "^•^^'^ 

Jn 



(15) 



^(plPLXSJT-jAjNJtAN) = ^[l-(-)'+^+^ 

X {plPLXSJJ\V\jAjNJtAN)- (16) 

T = {TMt.Mt = nij^ +mtj^}, Ian = {tA = ^/2,mt^ = -1/2, = l/2,mf^}, with 
nit^ = 1/2, and m;„ = —1/2, and / and L stand for quantum numbers of the relative 
and cm. orbital angular momenta in the AN system. The transition matrix elements 
^{plPLXSJT; JaJnJ^an) depend on the cm. and relative momenta, which are given 
in terms of the integration variable Ep by 



A-2. 



P = 2^ -^Mep, p = ^M{Aa^j, - ep), (17) 

where the energy conservation condition has been used. The angular momentum cou- 
plings 1+L = X, and A + S = J have been carried out, /= \/2J -\- 1, and A =Z+N +\ 
is the total number of baryons. 

It is self-evident that for A — oo one obtains the same result as in Refs. [6, 34, 36]. It 
is also worth noting that the overall outcome of the recoil on Fyy is very small, mostly 

because the effect of the factor (^^j^) ^"^^ in Eq. (14) is, to a great extent, cancelled by the 

effect of the factor {-^2) originating from yj ep(Aa^7^ — ep)dep. This is the reason 
why we have not included the recoil previously. 
From the relation 

2 

EnN = ^UnJf ~ 'T^Pi (18) 



which follows from (12) and (13), one can now easily derive the spectrum of Fat as a 
function of the sum energy E^^v [58]: 



4M^ 



M I /"^ ^ I 

^ ^0 rf.rTr. ^ 



(19) 



where 



I MA ( 



P = \Hr E-A[ 



ccnJf I ' 



p 



^2M{A-2){Aa^j,-E), (20) 

, A-2 
^unJf = ^«nJf ' (21) 



and the condition 



A'a^j,<E<Aa^j,, (22) 

has to be fulfilled for each contribution. 

In the same way from (14), and (15) we can easily arrive to an expression for as 
an integral on the cm. momentum P, namely 



with 



'' = yirr2sK^'' ''A.„.,=W^A„..„ (24) 



and similarly for P^. t is clear that the condition P < Pa^^j^^ has to be fulfilled for each 
contribution. 

Following step by step the developments done in Refs. [38, 58, 59], the Eq. (9) can be 
cast in the form 



= - XI [ dcosG f pljdpN I pldpn5{Aaffjp-ER-EnN)'^aNJF{pP)^('^^) 
^ UnJf J J J 



where the cm. and relative momenta, given in terms of the integration variables in (25) 
read 

P = ^{A-2){2MAa,JF-pl-pl), (26) 

and 



Next, the 5 function in (25) can be put in the form 
A -2 2M 



^[S{Pn-p^)+HPn-Pn)]. (28) 



Pn 



± - -(A-l)-i 



\Pn-Pn , 

where 

PN cos dnN T ^J2M{A-2){A-l ) A^^y^ - p2 [(A - 1 )2 _ cos2 

(29) 

The Eq. (9) becomes now 

= ITTTrr E /, dcos^nN PNdPN—r^\^a^JF{pP)p„-^pt 

+ (Pn^P;), (30) 

where the notation p„ — > indicates that ^{plPLXSJT; JaJnJ'^an) is to be computed 
with P and /? given by Eqs. (26) and (27) with pn replaced by p+. We have shown 
numerically that the last term in (33) is negligibly small in comparison with the first 
one, and therefore it will be omitted from no w on. 

With the simple change of variable p — > y/2ME one finally gets 



= (A-2)^£ r\co^enN l^'^'''' dEN\^En^a^j,{pP). (31) 
^ ^i„J-i Jo V E' 



a^Jp • 
where 

E' = {A-2){A-l)Aa^j^-EN[iA-l)^-cos^e„N], (32) 



En 



E' — WEn cos GnN 



(33) 



P = \l 2M{A - 2) (Aa^jp -En- En) , 



and 



p = ^ MAa^jp 



4(A-2) 



(34) 



(35) 



It might be worth noticing that, while E' does not have a direct physical meaning. En is 
the energy of the neutron that is the decay partner of the nucleon with energy En. The 
maximum energy of integration in (31) is 



pinax ^ ^ A 



a^Jp- 



(36) 



INDEPENDENT-PARTICLE SHELL-MODEL 



The formulas derived so far for the NMWD rates Fyy are totally general, and don't 
depend on the nuclear model that is used to describe the initial hypernuclear state |7/) 
and the final nuclear states | UmJf)- Beneath we describe the Independent-Particle Shell- 
Model (IPSM), which is widely used in finite nuclei. As usually done it will be assumed 
that the hyperon in the state j/^, with single-particle energy e^^, is weakly coupled 
to the A — 1 core, with spin Jc and energy Eq = Ejj — Cj^. Then the initial state is 
\Ji) = \ {JcjA)Ji), and 

Aa^jp=AM + Ec + £j^-Eaf,jp, with Ec = Ej, - (37) 
Moreover, the spectroscopic amplitude in Eq. (15) can be rewritten as 

{Ji\\{alal)j\aM = { } (7c||a)j|a^7^). (38) 
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FIGURE 7. Single particle energies of C. The particle states (above the Fermi level) lpi/2, and Id^ /2 
are the lowest states in '^N, while the hole states (below the Fermi level) lsi/2> and 1^3/2 are the lowest 
levels (in inverted order) in ^'B. 

Before proceeding a few words should be said on the s.p.e. of the emitted nucleon 
N = n,p in the two-nucleon decay An — )■ nN, which in the first section has been denoted 
by e^. Let us consider the decay j^^C— 7'^*^B+A + p, as one example. The s.p.e. of ^^C are 
displayed in Fig. 7. In the pure IPSM the particle states (above the Fermi level) lpi/2, 

and 1(^5/2 are the lowest states in ^^N, while the hole states (below the Fermi level) 
and 1 P3 /2 are the lowest levels (in inverted order) in ^ ^B, as illustrated in the upper panel 
in Fig. 8. The Is orbital is separated from the IP3/2 state by approximately 23 MeV, 
which is enough to break the 10 particle system of ^'^B, where the energy of the last 
excited state amounts to ~ 16.5 MeV. In fact, a single-particle state |ja?) that is deeply 
bound in the hypernucleus, after the NMWD can become a highly excited hole-state 
in the continuum of the residual nucleus. There it suddenly mixes up with more 
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FIGURE 8. Schematical representation of the l5i/2> and I/73/2 levels in ^'B within the models IPSM-a 
(upper panel), and IPSM-b (lower panel). Due to the coupling between hole and particle states, the deeply 
bound hole states acquire a width and is treated as continuum state in the IPSM-b model. 



complicated configurations (2hlp, 3h2p, . . . excitations, collective states, etc.) spreading 
its strength in a relatively wide energy interval [61], as schematically represented in the 
lower panel in in Fig. 8. ' Although the detailed structure and fragmentation of hole 
states are still not well known, the exclusive knockout reactions provide a wealth of 
information on the structure of single-nucleon states of nuclei. Excitation energies and 
widths of proton-hole states were systematically measured with quasifree {p, 2p) and 
(e, e'p) reactions, which revealed the existence of inner orbital shells in nuclei [62-70]. 

Therefore, the following two approaches for the final states \0CnJf) will be examined 
within the IPSM: 



IPSM-a 

Here, we completely ignore the residual interaction and, consequently, the only states 
\(Xj\/Jf) giving a nonzero result in Eq. (39) and therefore contributing to Eq. (16) are 
those obtained by the weak coupling, and properly anti symmetrizing, of the one hole 
(Ih) states to the core ground-state \ Jc)- Then, 

\ccnJf) ^ Ij'nJf) = \{JcJn^)Jf) , and £a^/^ ^ Ej^ = Ec- Cj^ , (39) 

where Cj^ is the single-particle energy of state j^, and the liberated energy in Eq. (37) 
becomes 



One should keep in mind that the mean life a A hyperon is Ta ~ 2.63 x 10 s, while the strong 
interaction times are of the order of 10^^' s. 



As an illustration, in the case of ^^Si the model space contains four single-particle 
states, both for protons and for neutrons (n^ = n„ = 4), namely, ^P3/2^ 

and 1^5/2- Thus, if the core state is \ Jc) — 11^^5/2"^^), the final states \ {JcJ]^^)Jf) are 
constructed by creating two holes in the ^^Si nucleus, and read: 

^^Si ^ nn + 26si 285^ ^np + ^^A\ 

|(lJ5/2n-i)2;0,2,4) |(lJ5/2n"ilrf5/2p-i);0, 1,2,3,4,5) 

|lJ5/2n"U5i/2«"i;2,3) \ld^i2n~'^lsij2P~^;2,3) 

1 1^5/2«~^ lPl/2«~^2, 3) I lJ5/2«~Hpi/2P~^;2, 3) 

1 1^^5/2«~ ^ lP3/2«~ ^ ; 1 , 2, 3 , 4) I IJ5 i2n- 1 1/73 /2P~ ^ ; 1 , 2, 3 , 4) . 

After making the substitution (40) in Eqs. (9)-(36) one can perform the summation on 
aj^Jp for each single-particle state jj^, as done in [6, Eqs. (11), (12), (13)], and do 

^aMpP)^^jN{pP) = E E \^{plPLXSJT;jAjNJtAN)\^, 



J=\jN-l/2\ SILXT 



(42) 



with the F^j are defined as 



Jf 

= ^-^Ej'; 'J^fmwjm^. (43) 



Jf 

The general formula (14), (19), (23), and (31), read now 
16M^ (A - 2\ 



jM^ ( A-1\ ' 1'^ I 

TiV = ^-^y-^) I d^pL\/£pi^h-£pWjNipP), (44) 

4M^ I ^ I 

rAr=^V^(A-2)3y^ dEnNY,^j{^jN-EnN){EnN-/S!j^)^j^{pP\ (45) 



2M I A -2 /-^A 

N = 



I 'dP%P^\IPl,-P^^mipP)^ (46) 

]N 



and 

Jn 

The meaning of all other quantities is self-evident from the initial expresions. 



Tn = (A-2)^J^ ^icos0„iv * dEN^^En^j^ipP). (47) 



IPSM-b 



Foraially, one starts from the unperturbed basis \iNJF)o with = 1,2, . . .nAr,niv + 
1, . . ., where for «V < nyv we have the same simple doorway states IJnJf) in Eq. (40) 
(listed in Eq. (41) for ^^Si), while for > t^n + ^ we have more complicated bound 
configurations (such as 3hlp, Ahlp, . . .in the case of ^^Si) as well as those including 
unbound single-particle states in the continuum. As in Ref. [61], the perturbed eigenkets 
I ocj^Jf) and eigenvalues Eaj^jp are obtained by diagonalizing the matrix o {InJf \H\i'j^JF)Q 
of the exact Hamiltonian H: 



with 



{ajsiJF \H\a'j^JF) = Eai^jp 5^^^'^^ (48) 



oo 

= LCff'\JNJF)+ I Cff'\iNJF)0. (49) 



matrix element (^c| |«/ \\0CnJf) in Eq. (10). Therefore, the Eq. (14) takes the form 



It is easy to see that only the ket IJnJf) in the expansion (49) will contribute to the 
16M3 /A-2^^/^ 

jnocnJf 



N 



I \ 3/2 „A I 



(50) 

where 



ILXSJT 



2 

• (51) 



To evaluate the amplitudes C^"^^ one would have to choose the appropriate Hamil- 
tonian H and the unperturbed basis \iNJF)o, and solve the eigenvalue problem (48). We 
will not do this here. Instead, we make a phenomenological estimate [58]. First, be- 
cause of the high density of states, we will convert the discrete energies A^^y^ into the 
continuous variable e, and the discrete sum on into an integral on e, i.e., 

/oo 
jCj,j,{e)\^pj,{e)de, (52) 

where py^(e) is the density of perturbed states with angular momentum In this way 
the Eq. (50) becomes 

Tat = 



fA — 2\' /"^ 

l^4~j J "^^^.^ /_^^/ePy^7^ (e)ep Vep(e - ep)^j^jp (pP), 

\ / ifjji7 °° 



JnJf ' 

(53) 



where 

Fj,j,ie) = \Cj,j,{£)\^pj,{e) (54) 

is called the strength function [61,71, 72] and represents the probability of finding the 
configuration IJnJf) = \ {Jc,jN^)'fF) per unit energy interval. Moreover, the Eq. (17) is 
substituted by 



P = 2y-^Mep, p=^/M{e-ep), (55) 

and the condition £ > £/> has to be fulfilled throughout the e integration. It is convenient 
to introduce the averaged strength function 

J Jc+Jn 

P^^(^) = HiwTTT ^ ^JnjA^)^ (56) 



where 



This allows to simplify Eq. (53) by making the approximation Fj^^jp (e)~Py^(e)toget 
16M3 /A-2\3/2 

JN 



N 



6M^ /A — 2 \ f°° 



(58) 

The IPSM-a results would be recovered if one made the further approximation 

¥j,{e) = d{e-^j,). (59) 

Here, in IPSM-b, the 5-functions (59) will be used for the strictly stationary states, while 
for the fragmented hole states we will use Breit-Wigner distributions, 

where are the widths of the resonance centroids at energies Ay^ (see [61, 
Eq.(2. 11.22)]). One proceeds similarly with the Eqs. (19), (23), and (31). It turns 
out that the expressions within the IPSM-b can be obtained from those of IPSM-a 
through the replacements: 

Aj^^e, and d£Fj^{e)---. (61) 

Jn Jn 



To evaluate the transition rates we need to know the spectroscopic factors F^j 

given by (43), which depend on the angular momenta Jc and Jj, experimental values 
of which are given in Table 1. It is also necessary to choose between the jj and LS 
couplings. As in the previous work [6] (see Table I) we used here the 77-coupling, 
which is extensively used in nuclear physics in view of large spin-orbit splitting. In fact, 
the experimental P1/2— P3/2 energy difference is ~ 6 MeV around ^^O. The resulting 
spectroscopic factors are shown in Table 2. 



TABLE 1. Experimental values of core spin Jc, and initial spin Jc for several A-hypernuclei. 


Nucleus 


iH iHe iH iHe iU ifie i'B i^C i^C i^N i^O i^O fSi 


1 1 


1 1/2 1/2 1 3 3/2 1 1/2 5/2 1 


1 1 


1/2 1 1/2 1/2 1/2 5/2 1 1/2 3/2 1 1/2 2 | 


TABLE 2. Spectroscopic factors F^j multiplied by 2j + 1. 


; NJ in iHe iH iHe iU ifie i^B i^C i^O i^O fSi 


si/2 nO 
nl 
pO 
/'I 


3/221111 1 111111 

1/2 033 3 3 3 33333 3 
3/2121111111111 

1/2 3 3 3 3 3 3 3 3 3 3 3 


P3/2 nl 
nl 

pi 


- - - 5/2 3 13/2 7 6 6 6 6 6 

- - - 3/2 5 11/2 5 10 10 10 10 10 

- - - 5/2 3 13/2 6 6 6 6 6 6 

- - - 3/2 5 11/2 10 10 10 10 10 10 


nl 
pO 
pl 


---------0011 

---------223 1 

---------0111 

---------233 3 


n2 
n3 
p2 
P3 


----- - - - -- --16 

______ _ _____i4 

______ _ _____i5 

______ _ _____21 



NMWD SPECTRA 

The transition probability densities Snwi^p), SnNi^nN) SnNiP), SnNi^os^nN)^ and 
Si^f{Ef^), can now be obtained by performing derivatives on Ep, EnN, P, Ej^j, and 
cos OnN in the appropriate equation for T^, namely, Eqs. (44),(45),(46), (47), and (47), 
respectively. 



Effects of the deeply bound hole states 

In Ref. [58] we have studied the effects of the deeply bound hole states on the 
correlation spectra SnN{EnN) of several hypemuclei. The Fig. 9 shows the normalized 
energy spectra Snp{E)/Tp for ^He, ^He, ]^C, ]^0, and ^^Si hypemuclei, evaluated 
within the full OMEP, that comprises the (;r, T] , ^, p , O), K*) mesons. Quite similar results 
are obtained for the nn pair, i.e., for S„„(i?)/r,j. The s.p.e.'s for the strictly stationary 
hole states have been taken from Wapstra and Gove's compilation [73], and those of the 
quasi-stationary ones have been estimated from the studies of the quasi-free scattering 
processes {p,2p) and {e,e'p) [62-70]. 
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FIGURE 9. (Color online) Normalized energy spectra Snp{E)/Tp for ^He, ^He, 
hypernuclei for the full OMEP obtained within the approaches IPSM-a (upper panel) and IPSM-b (lower 
panel). For the i-shell hypernuclei, only the IPSM-a approach has been used. 



^/^C, ^j^O, and ^^Si 



The two IPSM approaches exhibit some quite important differences: 

a) IPSM-a: The spectra cover the energy region 110 MeV < E < 170 MeV and con- 
tain one or more peaks, the number of which is equal to the number of shell-model 
orbitals l5'i/2, li?3/27 1/2 1 1 '^5/2 1/2, 1^^3/2 "■ ' that are either fully or partly oc- 
cupied in \ Jc). Before including the recoil, all these peaks would be just spikes at 
the liberated energies Aj^, as can be seen from (12) setting Er = 0. With the recoil 
effect, they behave as 

and develop rather narrow widths ~ [b^M{A — 2)]^^ where b is the harmonic 
oscillator size parameter, which has been taken from Ref. [40]. These widths go 
from = 3 MeV for ^^Si to = 20 MeV for ^He, as indicated in the upper panels of 
the just mentioned figures. 

b) IPSM-b: In the lower panels of the same figures are shown the results obtained when 
the recoil is convoluted with the Breit-Wigner distributions (60) for the strength 
functions of the fragmented deep hole states. The widths yj^ have been estimated 
from Refs. [61-70], and in particular from [62, Fig. 11] and [66, Table 1], with 



results: yi,,/^ = ^ MeV in ]^C, Yis^^^ = 14 MeV and 71^3^2 = 3 MeV in ]^0, ^ and 
Yisi/2 ~ 16 MeV and 71^3/2 = Yipi o ~ ^ MeV in ^^Si, both for protons and neutrons. 
One sees that, except for the ground states, the narrow peaks engendered by the 
recoil effect become now pretty wide bumps. 

We feel that the above rather rudimentary parameterization could be realistic enough 
for a qualitative discussion of the kinetic energy sum spectra. A more accurate model 
should be probably necessary for a full quantitative study and comparison with data. 



Interpretation of BNL experiment E788 on ^He 

Particularly interesting is the Brookhaven National Laboratory experiment E788 on 
^He, performed by Parker et al. [54], which highlighted that the effects of the Final 
State Interactions (FSI) on the one-nucleon induced decay, as well as the contributions 
of the two-nucleon induced decays, AA^A^ — > nNN, could be very small in this case, if 
any. 

Therefore one might hope that the IPSM could be an adequate framework to account 
for the NMWD spectra of this hypernucleus. This has been done in Ref. [60] by employ- 
ing the ;r + A" exchange potential, with soft cutoffs {K-n = 0.7 GeV and A^- = 0.9 GeV), 
which is capable of accounting for the experimental values related to Tp and r„ in all 
three ^H, ^He, and ^He hypernuclei [60]. It is labelled as SPKE model and is not very 
different from the PICE model used by Sasaki et al. [35]. 

The transition probability densities Sn{E), SnN{E), and S'„Ar(cos0) contain the same 
dynamics, but involve different phase-space kinematics for each case. In particular, the 
proton spectrum Sp{E) is related with the expected number of protons dNp{E) detected 
within the energy interval clE through the relation 



dNpjE) 
dE 



= Cp{E)Sp{E), (63) 



where Cp{E) depends on the proton experimental environment and includes all quanti- 
ties and effects not considered in Sp{E), such as the number of produced hypernuclei, 
the detection efficiency and acceptance, etc. In experiment E788, after correction for ac- 
ceptance, the remaining Cp{E) factor is approximately energy-independent in the region 
beyond the detection threshold, E^ [74]. In what follows, we will always compare our 
predictions with the experimental spectra that have been corrected for acceptance and 
take into account the detection threshold. Thus we can write, for the expected number 
of detected protons above this threshold, 

^P= I dE = Cp ' Sp{E)dE = CpTp. (64) 

Je^ dE "^Jeo 



The 3/2j peak is at 6.32 MeV, but small amounts of the P2/2 strength are also fragmented to the states 
of 9.93 MeV and 10.7 MeV [70]. 
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FIGURE 10. Comparison between the experimental and theoretical kinetic energy spectra for protons 
(upper panel) and neutrons (lower panel). The data are acceptance corrected [74], and the calculated 
results are obtained from Eqs. (67) and (68). 



This allows us to rewrite (63) in the form^ 

{E>El). (65) 



^Np(£) SJ^) 

dE r. 



p 



The spectrum Sp{E) is normalized to the experimental one by replacing Np in (65) with 
the acceptance-corrected number of actually observed protons, 



m 



N7 = £AN7(£,0, (66) 



where ANp^^(£',) is the acceptance-corrected number of protons measured at energy 
Ei within a fixed energy bin AE^, and m is the number of bins beyond the detection 
threshold. Thus, the quantity that we have to confront with data is 



mp{E)=Wp'P^p^^, (67) 



r 



p 



where the barred symbols (Nj/^ = 4546, and Tp = 0.168) indicate that the proton 
threshold E^ = 40 MeV [74] has been considered in the numerical evaluation of the 
corresponding quantities. In contrast to ANp'^^ (£■,), ANp(£') is a continuous function of 
E. 



A similar expression is valid for the j3-decay strength function (see, for instance, [75, Eq. (5)]). 




COS 

FIGURE 11. Comparison between experimental opening angle correlations for proton-neutron (upper 
panel) and neutron-neutron (lower panel) pairs. The data AN„^ (cos 0,) are acceptance corrected and do 
not contain events with En < 30 MeV [74]. The theoretical results are obtained from Eq. (70), with 
N"j^ only containing events with cos 9„n < —0.5. Two cases are presented: 1) Theory A, where both the 
angular and the single kinetic energy cuts are taken into account, and 2) Theory B, where the cuts are not 
considered in the calculations. 



As the one-proton (one-neutron) induced decay prompts the emission of an np (nn) 
pair, one has in the same way for the one-neutron spectrum 

AN„(£) = N^^'^AeJ-^^^^^^^^. (68) 

i p -|- zi „ 

Here, N^^^ = 3565, and Tp + 2r„ = 0.198 have been evaluated with a neutron thresh- 
old of 30 MeV [74]. In Fig. 10, our results are compared with the measurements of 
Parker et al. [54]. A similar, but somewhat different, procedure is followed for the co- 
incidence spectra. The main difference arises from the fact that the angular-correlation 
spectra, AN^^^(cos0/), as well as the kinetic energy sum data, AN^^(£',), besides be- 
ing acceptance-corrected, were measured with detection thresholds of 30 MeV for both 
neutrons and protons. More, in the selection of the kinetic energy sum data it was also 
applied an angular cut of cos 0„p/ < —0.5. In order to make the presentation simple, the 
observables that comprise only the energy cuts, and those that include both the energy 
and the angular cuts, will be indicated by putting, respectively, a tilde and a hat over the 
corresponding symbols. 

Thus, the number of nN pairs measured in coincidence can be expressed as 

n:^'' = lAN^Jlcose,) = £AN:J(£,), (69) 

/=1 i=l 



where the angular bins with cos di > —0.5 are excluded from the first summation. The 

- — exp — ~- exp 

AN„^ (cos 6j) and AN„^ (£■,) data should be compared, respectively, with 

AN„A,(cos0) =N:>cos0„^:H^, (70) 

and ^ 

^nNiE)=N:^^PAEj-^. (71) 

Here, from Ref. [74] N^^^ = 4821, = 2075, A cos G„n = 0.04 and AE„n = 10 MeV, 
while Tp = 0.1709 and r„ = 0.0113. These results (Theory A) are compared with the 
E788 data in Figs. 11 and 12. For completeness, in the same figures are also shown 
the results for 5„iv(cos 0) — )■ 5„Af(cos G), SnN{E) — )■ SnN{E) and Tn — )■ T^, i.e., when no 
energy and angular cuts are considered in the theoretical evaluation, and Tp = 0.1793 
and r„ = 0.0122 (Theory B). 
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FIGURE 12. Comparison between experimental kinetic energy sum spectra for proton-neutron (upper 
panel) and neutron-neutron (lower panel) pairs. The data AN„^ (£,) are acceptance corrected and only 
contain events with > 30 MeV and cos 0„n < —0.5 [74]. The theoretical results are obtained from Eq. 
(71), and two cases are shown: 1) Theory A, where both cuts are taken into account, and 2) Theory B, 
where the cuts are not considered in the calculations. 

We conclude that the overall agreement between the measurements of Parker et al. 
[54] and the present calculations is quite satisfactory, although we are not considering 
contributions coming from the two-body induced decay, ANN nNN, nor from the 
rescattering of the nucleons produced in the one-body induced decay, AN nN. How- 
ever, before ending the discussion we would like to point out that: 

1. As expected, the theoretical spectrum ANp(£'), shown in the upper panel of Fig. 
10, is peaked around 85 MeV, corresponding to the half of the 2- value Ap = 170 



MeV. Yet, as the single kinetic energy reaches rather abruptly its maximum value 
^mca _ ]^27 MeV (see Eq. (36)), the proton spectrum shape is not exactly that 
of a symmetric inverted bell. Something quite analogous happens in the case of 
neutrons, as can be seen in the lower panel of Fig. 10. The experimental data seem 
to behave in the same way. To some extent, this behavior of ANp(£') and AN„(£') 
is akin to the behavior of the AN„a?(£'), which suddenly collapse at the Q-values. 

2. There are no data at really low energies for the proton case which would allow to 
exclude the FSI effects for sure, and the neutron data for low energies are afflicted 
by large error bars. However, there is no need to invoke these effects, nor those of 
two-nucleon induced NMWD, to explain the data, as occurs in the proton spectrum 
of ^He [53]. This hints at a new puzzle in the NMWD, but it is difficult to discern 
whether it is of experimental or theoretical nature. 

3. The calculated spectra AN„p(cos0) shown in the upper panel of Fig. 11, are 

strongly peaked near 9 = 180", which agrees with data fairly well. However, while 
it is found experimentally that 28% of events occur at opening angles less than 120*^*, 
theoretically we get that only ^ 2% of events appear in this angular region. We find 
no explanation for this discrepancy. Nevertheless, the fact that not all events are 
concentrated at = 180", is not necessarily indicative of the contributions coming 
from the FSI or the AA^A^ — > nNN decay, as suggested in Ref. [54]. 

4. The calculated angular correlation AN„„(cos0), shown in the lower panel of Fig. 
11, is quite similar to that of the pn pair; that is, its back-to-back peak is very 
pronounced. This behavior is not exhibited by the experimental distribution. In 
addition, while 11% of events are found experimentally for cos0 > —0.5, in the 
calculation only ^ 3% of them appear at these angles. We feel however that, 
because of the poor statistics and large experimental errors, one should not attribute 
major importance to such disagreements. 

5. Both calculated kinetic energy sum distributions AN„a^(£'), shown in Fig. 12, 
present a bump at ?s 160 MeV, with a width of ^ 30 MeV, which for protons agrees 
fairly well with the experiment. We would like to stress once more that the spread- 
ing in strength here is totally normal even for a purely one-nucleon induced decay. 
The kink at ~ 130 MeV within the Theory A comes from the angular cut, and from 
this one can realize that the nN kinetic energy sum spectra below this energy are 
correlated with the angular coincidence spectra AN„Ar(cos0 < —0.5). The bump 
observed in the experimental AN„„(£') spectrum at 90 MeV is not reproduced 
by the theory, which may be indicative of nn coincidences originated from sources 
other than An decays, as already suggested in Ref. [17]. Another source for the 
difference between our model calculation and the data may be traced to np and 
nn final state interactions. Whereas in the former the intensity of this interaction 
is reduced owing to the Coulomb repulsion felt by the proton, in the latter the two 
neutrons may interact strongly and thus shift the peak to lower kinetic energy sum. 

In summary, to comprehend the recent measurements in ^He, we have outlined for 
the one-nucleon induced NMWD spectra a simple theoretical framework based on the 
IPSM. Once normalized to the transition rate, all the spectra are tailored basically by the 



kinematics of the corresponding phase space, depending very weakly on the dynamics 
governing the AN — t- nN transition proper. As a matter of fact, although not shown 
here, the normalized spectra calculated with the full PSVE model are, for all practical 
purposes, identical to those using the SPKE model, which we have amply discussed. In 
spite of the simplicity of the approach, a good agreement with data is obtained. This 
might indicate that, neither the FSI, nor the two-nucleon induced decay processes play 
a significant role in the i'-shell, at least not for ^He. 

OUTLOOK 

Before being detected the newborn nucleons in a NMWD suffer final state interactions 
(FSI) with the nuclear environment, and consequently their two-nucleon Fyv, and three- 
nucleon r2 decay rates are not observable from a quantum-mechanical point of view, 
as recently pointed out by Bauer and Garbarino [76]. These FSI give rise to emission 
of new secondary nucleons, that are counted by the detection systems together with the 
primordial ones, without being possible to distinguish ones from the others. The IPSM 
developed so far is a simple fully quantum-mechanical formalism for the theoretical 
investigation of decay rates Tf^ and their spectra. It don't describe neither the decay rate 
r2 nor the FSI. Therefore, it is not surprising that this model does not reproduce well the 
FINUDA experiment for the [53], as shown in Fig. 13, although it reproduces well 
[60] the BNL experiment for ^He [54]. 
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FIGURE 13. Comparison between the experimental [53] and theoretical [76] kinetic energy spectra for 
protons from ^^C decay. 

At the time being we are studying the proton kinetic energy spectra obtained in 
the FINUDA experiment for \ He, ^Li, ^Be, ]^B, ]^C, )^C, and [78, 79]. We 
are comparing their results with the simple IPSIVI, with the purpose to quantify the 
contributions of the FSI and the three-nucleon emission. The same is being done for 
the recent KEK measurments on angular correlations and kinetic energy sum of np and 
nn pairs [49, 50, 52, 80], as well as the cm. momentum spectra in [80]. Later on 
we will include the FSI, a consistent treatment of which would require in general a 
genuine three-body approach for the mutual interaction of the two emitted nucleons and 
the residual nucleus. Presumably due to the enormous computational challenges, this has 
never been tackled in the past. We are also planning to extend IPSIVI for the evaluation 



of the decay rate r2 for the emission of three primordial nucleons, which has been done 
so far only in the framework of the Fermi gas model [36, 42, 43, 81, 82]. 
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